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1 Introduction 



One of the most fascinating and technically demanding parts of the theory of two- 
dimensional integrable systems constitute the models with the spectral parameter on 
an elliptic curve. These include Landau-Lifshitz (LL) and Krichever-Novikov (KN) equa- 
tions, as well as elliptic Toda lattice (ETL) and elliptic Ruijsenaars-Toda lattice (ERTL). 
The latter two systems are less known; they will be discussed in detail in this paper (the 
word "elliptic" in their names refers to their relation to an elliptic curve). We will explain 
how all these models can be unified on the basis of a single equation introduced in [2] as 
the nonlinear superposition formula for Backlund transformations of the KN equation, 
and studied further in [HI EH]- We will denote it |Q4 following [Oj. This discrete model 
plays the role of the "master equation" in the (s/2 part of) 2D integrability: most of 
other models in this area can be obtained from this one by certain limiting procedures. 
Some of the interrelations have been known for a while, but several important parts of 
the picture were discovered only recently. The tower of continualizations is illustrated by 
the following diagram: 
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Here the left column denotes the independent variables domain of the system. The first 
line contains totally discrete systems on planar graphs studied recently in 013011113111, 
see also fSJEI]. The link between |Q 4 | and discrete ETL is a sort of reduction related 
with black- white coloring of quad-graphs, which was discovered in The transition 
r — ► Z 2 is just the specification of the general construction to the most important cases 
of regular lattices. Discrete equations in this more traditional setting were studied for a 
long time, see e.g. [H1I1IJII3I23I2I1I23I31II]- 

The arrows Z 2 — > Z x R — ► M? have a double meaning: a continuous system can be 
considered as a (non-classical) symmetry of the discrete one, or alternatively it can be 
obtained by a continuous limit. The amount of papers growths exponentially at each step 
down the diagram. We mention only the classical sources [2E1 H3 EI] and classification 
results on lattice systems due to the Siberian school j23 122 123 EE] • 

The contents of our paper correspond to the route along the arrows in the framed 
part of the above diagram (the rest part can be found, e.g. in [HIE]). In the main text we 
discuss only the non-degenerate case related to an elliptic curve. Several more simple cases 
can be obtained without additional efforts by some specification of the general formulae, 
for example by substituting sinh(x) or x instead of the Weierstrass function a(x). This 
corresponds to a degeneration of an elliptic curve into a rational one. All constructions 
of the present paper remain valid under such degenerations, which constitute the list 
IQ1HQ4I °f integrable discrete equations on quad-graphs placed in Appendix |XJ Actually, 
it contains six equations, since the parameter 5 in equations |Qi| and IQ3I can be scaled 
to either or 1. This list represents the most essential part of the classification result 
obtained in j^j. 



2 Discrete equations on quad-graphs 

We start with equations 

Q(u ,u 1 ,u 2 ,u 3 ;a,(3) = (1) 

on quad-graphs (that is, planar graphs with quadrilateral faces). Here the "fields" w fc G 
CP 1 are assigned cyclically to the vertices of any face of the graph and the parameters 
a,/? 6 C are assigned to its edges. It is required that opposite edges of any face carry one 
and the same parameter (see Fig. |TJ). According to [3 [THj, integrability of such equation 
is synonymous with 3D consistency, see below. The most important representative of 
integrable equations (JTJ) is the equation TJl which we give here in the nice form derived 



in [TB|: 

A((u - 6)(« 3 - b) - (a - b)(c - b)) ((«! - b)(u 2 -b)-(a-b){c- b)) 
+B((uo — a){u\ — a) — {b — a) (c — a)) ((7/3 — a) (1*2 — a) — (b — a)(c — a)) 

= ABC(a-b), (Q 4 ) 

where the points 

(a, A) = {p{a), p'(a)), (b, B) = p\(3)), (c, C) = (p(/3 - a), p'{(3 - a)) 

belong to the elliptic curve £ = {A 2 = r(a)}, and r(a) = 4a 3 — — g% is the notation 
for the Weierstrass polynomial used throughout the paper. So, in this case it is natural 
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Fig. 1: The basic quadrilateral and -D4 group of symmetry 



to assume that the parameters a, (3 belong to the period parallelogram of £ rather than 
to C. A collection of formulas related to 8 and to the corresponding Weierstrass elliptic 
functions is put in Appendix [Bl 

It is obvious that |Q4| is invariant under reflections in the diagonals: u <->■ w 2 , a ^ (3 
and u\ <-> U3, a «-> (3, as well as under reflections (w , W3) «-> (111,112) and {uq,u\) «-> 
(113, U2), so that this equation admits the symmetry group D4 of the square. 

Actually, the symmetry of is even more rich. This will follow from another form 



of this equation, which is of a fundamental importance also for many other reasons, as 
will be shown below. 



Proposition 1. Equation ~Ql is equivalent, under the point transformations Ui = p(xi), 
to equation 

F(x ,Xi;a)/F(x ,x 3 ;P) = F(x , x 2 ; a - (3), (2) 

where 

a(x + xi + a)a(x - x 1 + a) 

F(x 0j xi;a) = — ■ -. (3) 

a{XQ + x\ — a)a{xo — X\ — a) 

This proposition was proved in jH]. Equation (J2J is called the multiplicative three-leg 
form of (P) centered at Xq. Sometimes it is more convenient to use the additive three-leg 
form 

f(x Q , xi;a)- f(x , x 3 ; (3) = f(x Q , x 2 ; a- (3), (4) 

where ^ 

f(x Q ,Xi,a) = -logF(ar ,xi;a). (5) 

The functions / for equations IQiHQst as well as the corresponding point transformations 
of the fields and parameters are listed in Appendix [XJ 

Proposition ^ is checked with the help of some identities for elliptic functions listed 
in Appendix El With the help of (|B.4j) one sees immediately that the function is 
fractional-linear in m = p(xx) (but not in m = p(xo))\ 

( p(xi) - p(x + a) o- 2 (x + a) 

p(xi) - p(x - a) a 2 (x Q -a) 

Therefore, the formula (J2J yields an equation which is already amne-linear in ui,U2,u%, 
but with a complicated dependence on xq. Some additional transformations based on 
(jB.ljl . ()B.3j) are needed in order to show that this equation coincides with Q4 up to some 
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factor depending on xq. Note that eq. (J2J) defines, strictly speaking, a correspondence, 
since each Xi is defined up to a sign and up to shifts by periods of the function p(x). 
The function F has an obvious symmetry property 

F{Xo > Xl ' a) = F(zoX-,-<*) (7) 

which implies one of the -D4 symmetries X\ <-> £3, a <-> /?. However, another symmetry 
(xo,xs) <-> (xi,z 2 ) of the equation |Q 4 | becomes hidden in its three-leg form (J2J) centered 
at Xo- -Owe to this symmetry the three-leg form can be centered at an arbitrary vertex. 

On the other hand, three-leg form exhibits one more symmetry which remains hidden 
in the rational form IQ4I namely the fact that the diagonals of the quadrilateral are on 
the equal footing with its edges. Indeed, denote a = a% — 0.2, (3 = — Q2, then eq. (J2J) 
takes the form 

F(x , Xi] at - a 2 )F(x , x 2 ; a 2 - a 3 )F(x , x 3 ; a 3 - a x ) = 1. (8) 

This means that the full symmetry group of discrete KN equation is not the -D4 group of 
the square on Fig. Q but the symmetric group 5*4 of the tetrahedron on Fig. |21 Probably, 
there does not exist the form of equation, in which the whole group becomes apparent. 
Indeed, even in the simplest case of equation ( Q~i )s =0 (see Appendix [XJ) which reads 



CrossRatio(xo,xi,X2,xs) = a/ (3, this full group S4 is encoded in the transformation 
properties of the cross-ratio. 

As the first application of the three-leg form, we will prove that it (together with the 
above mentioned symmetry properties) implies the fundamental property of 3D consis- 
tency, which can be accepted as the definition of integrability for equations of the form 
((H). Recall [3 EE] that equation (JTJ) is called 3D-consistent, if the following holds. For 
arbitrary initial data Woj u i, U2, W3, define the values uu, U23, W31 (the enumeration, which 
is different from the one used elsewhere in this paper, is explained on Fig. EJ) according 
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to the equations 

Q(u ,u 1 ,u 12 ,u 2 ;a 1 ,a 2 ) = 0, (9) 

Q(u ,U2 } U23,u 3 ;a2,a 3 ) = 0, (10) 

Q(u ,u 3 ,u 31 ,ui;a 3 ,ai) = 0, (11) 

which correspond to three faces of the cube adjacent to the vertex 0. Then equations 

Q(u 1 ,u 12 ,u 1 23,u 31 ;a 2 ,a 3 ) = 0, (12) 
Q(u 2: u 23: Ui 23 ,u 12 ;a 3 ,ax) = 0, (13) 
Q(u 3 ,u 31 ,ui2 3 ,u 23 ;ai,a 2 ) = 0, (14) 

corresponding to the rest faces, define one and the same value U\ 23 . 

Proposition 2. Suppose that equation (OP admits the D± symmetry group and is equiva- 
lent, under some point transformation u = 4>{x), a = p(a), to the three-leg equation (0|). 
Then it is 3D- consistent. 

Proof. It is enough to show that if u\ 23 is defined by equation (fT2"|) . then (fT3~|) is fulfilled as 
well. Rewrite equations (jHJ), (JUJ), (fT2"|) (corresponding to the faces adjacent to the vertex 
1) in the three- leg forms centered at X\ (recall that this is possible due to the symmetry 
properties): 

F(x 1 ,x 31 ;a 3 )/F(x 1 ,x ;a 1 ) = F(x 1 ,x 3 \a 3 - ot x ), 
F(xi, xi2] a 2 ) I 'F(xx, x ; «i) = F(x 1 ,x 2 ;a 2 - ax), 
F(x 1 ,xi 2 ; a 2 )/F(x 1 ,x 31 ; a 3 ) = F(x 1: x 123 ; a 2 - a 3 ). 

From these there follows the equation relating the fields at the vertices of the dashed 
tetrahedron on Fig. EJ 

F(xi, x 2 ; a 2 - a x )/F{x u x 3 ] a 3 - a x ) = F(xt, x 123 ; a 2 - a 3 ). 
This is nothing but the three-leg form of the equation 

Q(ui, u 2 , u 3 , u 123 ; a 2 - a 1} a 2 - a 3 ) = 0, (15) 
centered at x\. This can be centered at x 2 , as well, resulting in the cyclic shift of indices: 

F(x 2 , x 3 ; a 3 - a 2 )/F(x 2 , xi, a x - a 2 ) = F(x 2 , x 123 ; a 3 - a x ). 

The latter equation, together with the three-leg forms of equations (jHl), (fTUj) centered at 
X2, leads to the three-leg form of (fT??|) . as required. □ 

3 Discrete Toda systems on graphs 

Existence of the three-leg form allows us to establish a direct link to discrete Toda systems 
on graphs introduced in 0. This link was discovered in j2j and is described as follows. 
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Fig. 4: Faces adjacent to the vertex xq. Fig. 5: The black graph star of the ver- 
tex Xq. 



Suppose that the quad-graph T is bi-partite, i.e. its vertices can be colored black and 
white so that the ends of any edge are of different colors. Consider the graph To whose 
set of vertices ^(r ) consists of the black vertices of T, and whose edges are diagonals 
of quadrilateral faces of T. Consider a vertex xq G ^(r ) (for brevity we will identify 
vertices with the field variables sitting there), and let Xq be a common vertex of n adjacent 
quadrilateral faces (xo, x 2 j-i, X2j, X23+1) G F(T), j = 1, . . . , n, like shown on Fig.HJ Then, 
multiplying the equations of the type (j2J) centered in x for all these faces, we see that 
the white vertices £2j-i cancel out and the black vertices X2j satisfy the equation on the 
star of the vertex Xq: 



Y[F(x ,x 2j ;a j 



— a 



j+V 



(16) 



or, in the additive form, 



(17) 



This is a discrete Toda system on the graph r . (Of course, a similar Toda type system 
holds also for the white graph Ti which is dual to r .) The parameters aj which were 
assigned originally to the edges (xo,X2j-i) G E(T), label now the corners of the faces 
of r . Notice that only the differences of these parameters appear in equations (fTTj) . so 
that simultaneous shifting all the parameters of the equations on the original quad-graph 
a — > a + A leads to the same Toda system on r . This shift plays the role of the spectral 
parameter in the zero curvature representation of the Toda system. 

As established in (3 QHJ, any 3D consistent equation (jlj is its own zero curvature 
representation. We will show that it plays this role not only for itself, but for a whole 
set of related systems, including discrete Toda systems. For this aim, rewrite eq. on 
the quadrilateral (u , ui, u 2 , u 3 ) G F(T) as 

k U UiU 2 U3 + ki(uiU 2 U 3 + UoU 2 U 3 + U UiU 3 + U UiU 2 ) + k 2 (u U 2 + U1U3) 

+ k 3 (u ux + u 2 u 3 ) + k A (u u 3 + uxu 2 ) + k 5 (u + u l + u 2 + u 3 ) + k 6 = 0. 



6 



(Concrete expressions of the coefficients ki in terms of (a, A) and (b, B) can be found in 
the original paper j2J.) The field 1*3 is a fractional-linear function (Mobius transformation) 
of U\ with coefficients depending on m , «2 and a,\, a 2 . Denote the matrix of this Mobius 
transformation as 

, _ . _ J_ f-k 1 u u 2 - k 3 u - k A u 2 - k 5 -k 2 u u 2 - k 5 (u + u 2 ) - fc e 

M 2, on, at 2 ) p y k u u 2 + ki(u + W2) + h kxu u 2 + k 4 u + k 3 u 2 + k 5 

(18) 

where D is the normalizing factor such that detC/ = 1. Considering now n quadri- 
lateral faces (u , u 2 j-i, u 2 j, ^2.7+1) £ F(r) adjacent to w , we see that, according to the 
cancellation property leading to the Toda systems, the product 

\ I U(u ,u 2j ;aj,a j+1 ) 

l<i<n 

generates an identical Mobius transformation and is therefore a scalar matrix. Moreover, 
the latter property is equivalent to the discrete Toda equation at the vertex u G V(T ), 
formulated now in terms of the variables (uq, u 2 j). Thus, we see that the following holds. 

Proposition 3. Discrete Toda equation Mb}) at the vertex uq G V^(ro) is rational in 
terms of uq and u 2 j, j = 1, . . . ,n. The discrete Toda system on an arbitrary graph To 
admits a zero curvature representation in the sense of J^, with the transition matrices 
across the edge (u , u 2 ) G E(T ) given by U (u , u 2 ; a\ + A, a 2 + A) . 

The fact that the Toda equation for the star of the vertex m is rational with respect 
to u 2 j = p(x 2 j), j — 1, . . . , n, is easy to see from (fTK|) . which can be rewritten with the 
help of the formula (J0J) as 

t~t p(x 2 j) - p(x + otj - a j+ i) a 2 (x + a>j - a j+1 ) _ 
~ = i p( x 2j) ~ p(xo - otj + a j+ i) cr 2 (x - aj + a j+1 ) 

However, the rationality with respect to Uq = p(xq) remains hidden in this representation. 
This is a new manifestation of the fact that the original (polynomial in Uj) and the three- 
leg forms of the equation Q4 are complementary in showing and hiding various properties. 



The next, Lagrangian property is better visible again in the variables x, with the help 
of the three-leg form. It is easy to see that the function df{xo,xi,a)/dxi is symmetric 
with respect to the flip xq <-> x±. Integrating this function with respect to Xq and x±, we 
see that there exists a symmetric function A(xq,xi; a) = A(xi,xo; a), such that 

f(x , xx\ a) = dA(x , xi, a)/dx . 

Proposition 4. Let E(Vq) be the set of the edges of the black graph Tq. Let the pairs 
of parameters be assigned to the edges from E(T ) according to Fig. [3J so that, e.g., the 
pair (ai,a 2 ) corresponds to the edge (x ,x 2 ). Then the discrete Toda equations fT7| ) are 
Euler-Lagrange equations for the action functional 

S— ^2 A(x ,x 2 ; en - a 2 )- 

{x ,x 2 )£E(r ) 
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Finally, we give the matrix formulation for the cancellation phenomenon which is the 
main feature of the three-leg form. 

Lemma 5. If the relation 

f(x ,xi;a) - f(x o ,x 3 ;0) = q <^> F(x , X\\ a)/F(x Q , x 3 ; 0) = e 2g (19) 

holds for some quantity q, then u 3 = 80(2:3) is a Mobius transformation of u\ = p(x±) 
with the matrix 

L = e q N(x ; 0, a) - e~ q N(x ] -0, -a), (20) 

where 

= a\x + 0)a\x - a) f-p(x + 0) p(x + 0)p(x - a)\ 

v 11 ' ' a(2x) \ - 1 p(x-a) J K J 

Proof. According to (JBJ), the formula is equivalent to: 

p(x 3 ) - p(x + (3) a 2 (x o + 0) p(xi) - p(x Q + a) a 2 (x + a) 



p(x 3 ) - p(x - 0) a 2 (x o -0) p(xi) - p(x - a) a 2 (x - a) 



Now it is a matter of a straightforward computation to represent this as the Mobius 
transformation from p(x\) to p(x 3 ) with the matrix L given by (12 OH . (jS} - □ 

Of course, matrix L is defined only up to a scalar factor, which is chosen in (|2U|). (|21|) 
so that detL = a (2a) a (20) = const. The telescopic cancellation taking place by adding 
several equations of the type (|19jl. is translated for the matrices L into the following 
property: if, additionally to (pUj). we have 

L x = e^N(x 0] 7, 0) - e-^N(x 0] -7, -0), 

then 

L\L = a(20) (e^ +q N(x ; 7, a) - e^A^o; -7, -a)) • 
This follows easily from the properties of the rank 1 matrices N, namely: 

N(x; 7, 0)N(x; 0, a) = a(20)N(x; 7, a), N(x; 7, 0)N(x; -0, -a) = 0. 

Notice that the case q = f(x ,x 2 ;a — 0) in eq. (fT9*j) corresponds to the additive 
three-leg form (0J) of equation |Q 4 | In this case we obtain from Lemma El a new formula 
for the transition matrices from Proposition EJ 

U(uo, u 2 - a x , a 2 ) ~ e f{x °' X2m ~ a2) N(x ; a 2 , ax) - e - f{xo ' X2m ~ a2) N(x ; -a 2 , - ai ), 

where symbol ~ means equality up to a constant scalar factor. A direct verification of 
this formula is actually equivalent to a proof of Proposition ^ 
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4 Lagrangian systems with discrete time 
and space-time slicing of regular lattices 

Now we consider the specification of the general scheme corresponding to the case of the 
shift-invariant graphs, or lattices. Recall the general construction of the discrete time 
Lagrangian mechanics JH]. Let X be a vector space (the theory can be developed also 
in a more general setting, when X is a manifold, but we will not need this here). Let 
£ : X x X ~R be a, smooth function, called the discrete time Lagrange function. For 
an arbitrary sequence X : eL D {ae, as + £,... , be} \— > X one considers the discrete time 
action functional 

6-1 

S = £(X{ne),X(ne + e)). 

n=a 

The sequences serving as critical points of the action functional (in the class of variations 
preserving X(ae) and X(be)) satisfy the discrete time Euler- Lagrange equation, which we 
shall write in the index-free form as 

dC(X,X)/dX + 8C(X,X)/dX = (22) 

where X, X, X stand for X(ne — e), X(ne), X(ne + e), respectively. This is an implicit 
equation for X. In general, it has more than one solution, and therefore defines a cor- 
respondence (multi- valued map) (X, X) i— > (X, X). To discuss symplectic properties of 
this correspondence, define 

P = OC(X,X)/dX e T* X X. 

Then (I22j) may be rewritten as the system 

P = -dC(X, X)/dX, P = dC(X, X)/dX. (23) 

This system will be called the Hamiltonian form of eq. (|22jl. It defines a multi- valued 
map (X, P) (X,P) on T*X. More precisely, the first equation in f!23l) is an implicit 
equation for X, while the second one allows one to calculate P explicitly and uniquely, 
once X and X are known. The fundamental property of this multi-valued map is the 
following: each branch of the map T*X \— > T*X defined by H3J) is symplectic with respect 
to the standard symplectic structure onT*X . 

The continuous limit of this construction is as follows. Suppose that X(ne) approxi- 
mates the smooth function X(t) at t = ne with some small e > 0. Suppose that by the 
substitution (X, X) = (X, X + eX + o(e)) there holds an asymptotic relation: 

£(X,X) ^ £ L(X,X) 

with some smooth function L : TX \— > R, where the sign pa means equality up to additive 
terms o(e). Then we have: 

OL/dX pa dC(X 1 X)/dX = P, 

OL/dX Pa e- 1 (dC{X,X)/dX + dC{X,X)/dX) = e-\P - P). 



9 




X X\ 

Fig. 6: Slice (4 + 2) on triangular lattice 



Thus, the discrete Lagrangian equations (|23|) approximate (serve as a discretization of) 
the continuous Euler-Lagrange equations with the Lagrangian L(X,X): 

P = <9L(X, X)/dX, P = dL(X, X) jdX. 

Now we can introduce the discrete time Lagrangian formalism for discrete Toda sys- 
tems on some regular two-dimensional lattices. We will consider planar graphs T with 
the set of vertices V(r) = Z 2 only. They will be distinguished by their sets of edges 
E(T). It will be convenient to think of the vertices (k, n) G Z 2 as representing the points 
(k,ne) G M 2 with some (small) e > 0. The first index k will enumerate the lattice sites, 
while the second index n has the meaning of the discrete time t = ne. For a fixed t = ne, 
we will write V 1 = {(k,t) : k G Z}, so that V(T) = \JteeZ Vt - 

A slicing of a graph T consists in choosing subgraphs T* so that V(r*) = V l ~ e U V f , 
the sets S(r*) = E l are disjoint, and E(T) = [j^E 1 . 

The configuration space of our discrete time systems is X = {X : Z i— > M.} consisting 
of real-valued sequences X = (xk)k&- An element P G T^X is a sequence P = {pk)k£Z- 
Clearly, for any fixed t = ne, the set of functions on V 1 is identified with X. Functions on 
V(r*) are naturally identified with pairs of functions (X, X) on V l ~ e and V*, respectively, 
and therefore belong to X x X. The discrete time Lagrange functions for our Toda systems 
are of the form 

C(X,X)= A(x,x';a-/5), 

i.e. are obtained by summing up the elementary Lagrangians over the edges of one slice 
of the graph T. 

The pictures of the slicing for the graphs playing the main roles in our subsequent 
presentation are given below. Starting from this point, we use the following abbreviated 
notation for the sequences like (xk)kez : we write x, x±i for Xk, Xk±i, respectively. This 
notation has not to be confused with indices used in the previous section. 

The central role in what follows will play the discrete Toda system on the triangular 
lattice and on the square lattice. All vertices of the regular triangular lattice have valence 
6, so that the corresponding Toda system is a 7-point scheme. Similarly, the Toda system 
on the regular square lattice is a 5-point scheme. The Hamiltonian properties of these 
systems depend essentially on the slicing chosen. For the triangular lattice we will choose 
the slicing on Fig. El The corresponding map (X,X) i— ► (X, X) is defined implicitly, 
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since each equation contains two fields x and x-\ from the time level t = ne + e). For the 
regular square lattice we will use two different slicings. The first one is shown on Fig. [7| 
in a somewhat unusual way (by a standard drawing the slice would run diagonally). It 
corresponds to an implicit discrete Toda system. Indeed, also in this case each equation 
contains two fields x and x_i from the time level t = ne + e. Finally, the slicing shown on 
Fig. |H1 represents an explicit discrete Toda system, since each equation contains only one 
updated field x. These two types of 5-point equations can be obtained from the 7-point 
one by reduction consisting in erasing some edges of the triangular lattice. 




X Xi x 



Fig. 7: Slice (2 + 2) on square lattice Fig. 8: Slice (3 + 1) on square lattice 



5 Triangular lattice 

Time discretizations of lattices of the Ruijsenaars-Toda type j2l] were introduced in 
[231 Ell ■ They were identified in [3] as discrete systems on the regular triangular lattice. 

The regular triangular lattice T can be considered as the black subgraph of the quad- 
graph known as the dual kagome lattice (drawn on Fig.|§]in dashed lines). The latter graph 
has vertices of two kinds, black vertices of valence 6 and white vertices of valence 3, and 
edges of three types, all edges of each type being parallel. The parameters corresponding 
to these three types of edges are denoted 0:1,0:2, 0:3 as shown on Fig. El It is easy to 
see that the parameters 03 are constant along horizontal stripes, that is 0-3 = a^, n and 
analogously, a 2 = o 2 ,fc and a,\ = oii >k+n . 

In what follows, we restrict ourselves to the particular case which admits a well-defined 
continuous limit, when the parameters ot\ and a 3 are constant: 

«i = A, a 2 , k = A - p k , «3 = A - e, (24) 

where A will play the role of the spectral parameter in the zero curvature representation of 
the Toda system. The discrete Lagrange function obtained by summing up the elementary 
Lagrangians along all edges of a time slice, is equal to 

C(X,X) = y^ y (A(x k , x k \e) + A(x k , x k+1 \(3 k - e) + A(x k , x k+1 ; -(3 k )). (25) 

k 

Equation of the discrete elliptic Toda system on the triangular lattice reads: 
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X 




Fig. 9: Fields and wave functions on the triangular lattice 



f(x, xi, -ft) + f(x, x; e) + f(x, z_i; /3_i - e) 

+ f(x, x-i', -ft-i) + f(x, x; e) + f(x, x 1 ;ft-e) = 0. (26) 

Using the fact that f{x,x';—a) = —f(x,x';a), we put the Hamiltonian form of this 
equation as: 

V = -f(x, x; e) + f(x, xi, ft) + f(x, x_i; ft-i) - f(x, X- X ; /3_i - e), (27) 
p = f(x,x;e) + f(x,x 1 ;ft-e). (28) 

Proposition 6. The system \2T\ l, h2ty) admits a zero curvature representation 

LV = V\L (29) 

with the matrices 

L = L(x,p, ft; A) = e p M(x, ft; A) + e~ p M(x, -ft; -A), (30) 
where the matrix M(x, ft; A) = N(x; A — ft, A)/cx(2A) zs gwen by 

v ' M ' ; <r(2a;)<7(2A) \ _1 - A) J 

The matrices V in \29\) are given by 

V = V(x, ar_i, e; A) = e eG M(x, e; A) + e~ eG M(x, -e; -A), (32) 

where 

eG = f(x, x_i; - /(x, /3_i - e). (33) 
12 



Proof. Consider the full system of the equations (JTJ) on the dual kagome lattice and denote 
the vertices of the white sublattice as shown on the Fig. El The auxiliary linear problem 
= corresponds to crossing two edges (x, x) and (x, X\). Adding the three-leg 
equations on the quadrilaterals (x, \l>, x, $1) and (x, $1, x%, \I/i), and using the formula 
(|2E1) we find: 

f(x, vl/; A) - f(x, * i; \-P) = f{x, x- e) + /(x, x 1 ;/3-e)=p. (34) 

Similarly, the auxiliary linear problem \1> = corresponds to crossing two edg 

and [XjX-i). To find V, add the three-leg equations on two quadrilaterals (x, <E>,x_i, ty) 

and (x, x_ 1; $): 

/(x,tt;A)-/(x,tt;A-e)=eG, (35) 

with eG defined by (|3*3*|). Now Lemma applied to (|3*4*j) . (|3*5|) . yields the formulas (|3*0|l. 
respectively □ 

Of course, in the Hamiltonian framework the main role in the zero curvature repre- 
sentation belongs to L, the transition matrix along the lattice. For instance, the so-called 
monodromy matrix Y\ k Lixk,Pk, Pk] A) is a generating function of integrals of motion. In 
this context it is important that Lk is local, i.e. depends on one pair of canonically con- 
jugate variables Xk,Pk only. As for the matrix V, it specifies just one (discrete) flow of 
the hierarchy. In the next section we will present another (continuous) commuting flow. 



6 Elliptic Ruijsenaars— To da lattice 

The Lagrangian (|2l)j) has a well defined limit at e — > 0, in the sense that £(X, X + eX) « 
eLiX,X). Alternatively, this limit can be performed directly in the equations of motion 
(J2ZJ), (EH) - Taking into account the relations 

_ x + 1 1 
lim fix, x; e) = log -, lim x; e) + /(x, x; e)) = 2((2x), 

e^O X — 1 e->0 e 

one finds, as the limit of the equation (|25|) and the one obtained by subtracting (|27|) from 



x + 1 

P = log- + fix,xx;P), (36) 

x — 1 

.dfix,xr,p) . dfix,X-i,P-i) <9/(x, x x ; /3) <9/(x, x_ x ; /3_ a ) 

p = x & + x ~ l dxZ, dp dju + C( x) - 

(37) 

From here the Newtonian equations of the elliptic Ruijsenaars-Toda lattice follow: 

x . dfix,x x ;P) . df(x,x-i]P-t) 
- X\ x_ r 



x 2 — 1 dxi 9x_i 



dfjx,xi,P) dfjx^.uP^) 

dp + WIi zazxh [6t 
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where 

?l^liS. = ^{x + x 1 + P)-C(x-x 1 + P)-C{x + x 1 -0) + C(x-x 1 -fi)), 

^01^1 = i(c(x + x x + (3) + C(x - x x + (3) + C{x + x x - (3) + C(x - x x - (3)). 

A remarkable feature of the zero curvature representation of Proposition El is that the 
matrix L in (|3U|) does not depend on e. Therefore the flow (|36|h (|37p shares the matrix L 
with the map (|27j) . ()28|) . In other words, the latter is a Backlund transformation of the 
former. 

Proposition 7. Elliptic RTL l&b}) . fiTfy admits a zero curvature representation 

L = A\L — LA (39) 

un£/i L from fifty and A = A(x,X-i,X-i\ A) = dV/de\ £=0 , where V is given in / f3^)) . 

By the usual substitution u = p(x), the elliptic RTL (|3*8*j) is brought into a rational 
form (which appeared first in [I]): 

2u — r'(u) iii u-i d . . m _. n w / in \ 

at^t = " M^-^y + k(u u \ b a + Yu log tti; tt - i; ' (40) 

where h is a biquadratic polynomial with discriminant r(u) defined in Appendix El By 
the derivation of eq. ()40|) from (|38|) the following alternative expression is used: 

2 df(x, Xl ;P) a(2/3)a(2x)a(2x 1 ) 

dxi a(x + Xi + (3)a{x — X\ + (3)a(x + x x — (3)a{x — X\ — (3) 

= p\x)p\x l ) 
h(u,ui,(3) ' 

It is obtained with the help of (|B.2jl . ()B.6j) . (|B.7jl . As a consequence, one finds: 



2 df^x_±P) = d dfjxx^d) + 
op OX OX i 

= p'(x)-^logh(u, m- (3) - + 2((2x), 

which is the second expression used in the derivation of eq. (|40p. 

7 Triangular lattice in Z 3 

Proposition 8. Consider a solution x : V(T) i— ► C of the discrete Toda system \2b}) on 
the triangular lattice T . Define a function y : V(T) i— > C by equations 

F{x, x- e) = F{x, y- l3)/F{x, x x -(3-e). (41) 
Then also the following equations hold: 

F(y } y; e) = F(y } x; (3)/F(y, y_ t - (3 - e). (42) 
Moreover, y is also a solution of the discrete Toda system Wty . 
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Fig. 10: Map x \— > y 



The geometric description of equation ()41)1 is the following: the variables y should be 
thought of as defined on the vertices of the new copy of triangular lattice which correspond 
to the half of the faces of the original one, namely those shaded on the Fig. El (so, y sits 
on vertex dual to the triangle (x, x, Xi)). Then equation (jl^j) reflects the fact that the 
variables x are attached to the half of the faces of the ^/-lattice. Moreover, equations (@H), 
f)42jl are nothing but the three-leg forms (JEJ) of the basic equation for the tetrahedra 



(x, x, xi, y) and (y,y,y-i,x). In terms of u = p(x), v = p(y), the equations (jUJ), (J32J 
read: 

Q(u, u, Mi, v] e, /3) = 0, Q(v , v, v_i,u; e, (3) = 0. 

Proof. Consider the system on the cubic lattice Z 3 = {(11^2,13)} consisting of equation 
|Q4| on each square face, with the parameters a m attached to the edges of the direction 
i m . This is possible due to the 3D consistency of equation |Q 4 | established in Proposition 
El Vertices of the x- and ^-lattices can be interpreted as the points of the cubic lattice 
lying in the planes ^1 + ^2 + ^3 = and ^1 + ^2 + ^3 = 2, respectively, see Fig.EJ Due to the 
three- leg property, both fields x and y satisfy the discrete Toda system (j2IIJ). Equations 
dHJ), are nothing but the tetrahedron equations (fT5]l for the corresponding cubes. 

Conversely, it is easy to see that, given a solution x of (JH3), one reconstructs uniquely 
the black sublattice of the cubic lattice, and then an arbitrary value at one white vertex 
determines the white sublattice uniquely. □ 

So, equations (jUJ), (|4~2|) can be considered as a sort of Backlund transformation x 1— > y 
of the discrete Toda system (pj)) . On the other hand, they can be considered as a 
discrete time system (mapping) (X, Y) i-> (X,Y), where X = (x k ) k& , Y = {y k ) k&L . In 
this interpretation, this mapping serves as discretization of the Shabat-Yamilov lattice 
considered in the next section. 

Proposition 9. The discrete time system \41\ ), \4ty admits a zero curvature represen- 
tation \2ty) with the matrices 

L{x,y,P;\) = U(u,v;\,\-P), (43) 
V(x,y„ u e;\) = U(u, £L i; A, A - e), (44) 

where u = p(x), v = p(y). 
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X 



Fig. 11: Embedding of the triangular lattice into the cubic one 



Proof. The matrices L, V correspond, as before, to the transitions ^ — > and ^ — > 
ty. As can be seen from Fig. El these are transitions across diagonals of square faces 
(x, y, ^i) and (x, y-i, of the cubic lattice. □ 

8 Shabat-Yamilov lattice 

The continuous limit e — > in equations (}4*T]) . (|42|) leads to the differential-difference 
equations 

x + 1 _ F(x,y;l3) y - 1 _ F(y,x;(3) 

x-1 F{x,x x -py y + 1 F(y,y- X -P)' 1 } 

It should be noticed that elimination of y from these equations yields ERTL ([38)1 . 

Under the change of variables w = v = p(y), the system (J45|) is brought into a 

rational form: 

2p(u 1 ,u,v;/3) . 2p(u,v,V-i,/3) 
u = , v = (46) 

«1 — V U — t>_! 

where the polynomial p is given in Appendix [Bj in the derivation of ()46|) from (|4"5)l one 
uses the formula 

F(x,y;/3) = p(y) - p(x + (3) _ p{x l ) - p(x - (3) 
F(x,xi,(3) p(y)-p(x-[3) p(x x ) - p(x + /?) ' 

and then takes (jB.8|) into account. 

The system (j4*HJ) was discovered in [221- Its relation to the semi-discretization of the 
LL equation is discussed in 

Proposition 10. The system \4b\) admits a zero curvature representation \3fy) with the 
matrices L from UJfy and A = A(u, u_i; A) = dU(u, t>_i; A, A — s)/de\ £= Q. 
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9 Elliptic Volterra lattice 



The same system on the cubic lattice that produced in Sect. 0the discrete time Shabat- 
Yamilov lattice serves as the origin of another important lattice system and its time 
discretization. To show this, consider one row of elementary cubes and denote the fields 
variables as on Fig. The parameters a, take the same values (|24|) as before. Note 
that in these notations a 2 (and therefore (3) are site- independent. Then it is easy to see 
that the fields x and x satisfy the equation 



F(x, x; e) = F(x, x\\ (3)/F(x, X-\\ (3 — e) 



(47) 



This equation is nothing but the three-leg form (jHJ of the basic equation for the 
tetrahedron {x, x±, x, It defines an (implicit) map X t— > X, where, as usual, X = 

x-i \j/ xi x 3 



«3 = A — e 




& 



.. ■■••''a 2 A > 



x 2 «i — A x 4 

Fig. 12: Discrete time Volterra lattice 



Proposition 11. The map OTp admits a zero curvature representation \29§ with the 
matrices 

L = L(xi, x, (3; A) = U(ui, u; A, A — (3), 
V = V(xi,x-i,e; A) = U(ui, A, A — e). 

Proof. As usual, the matrix L corresponds to the transition ^ — > ^1, the matrix V 
corresponds to the transition ^ — >• One sees readily from Fig. El that these are 
transitions across diagonals of elementary squares. □ 

In the limit e — > the equation ()47|) turns into the elliptic Volterra lattice 

x - 1 = F(x,Xi,P) 

x + l F(x,X-i;P)' 1 ' 

In the rational form, 

2p(u 1 ,u,u_ 1 ;p) 

u = , 

U\ — W_i 

it is due to Yamilov [22] . This lattice admits a zero curvature representation (}39|) with the 
same matrix L as in the discrete time model. In particular, the map (|4"7j) is a Backlund 
transformation for the flow (|48jl . 
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10 Square lattice, implicit scheme 



We consider now the skew square lattice with a slicing as on Fig. obtained from the 
triangular one by erasing the horizontal edges. This can be achieved by setting a± — «2, 
so that for our particular setting (|24jl all parameters (5 k of the Toda lattice vanish, and 
the parameters of the corresponding quad-graph (which is a finer square lattice shown 
by dashes on Fig. [T3J) are a x = A, a 3 = A — e. 




Fig. 13: Fields and wave functions on a skew square lattice 

All results on the discrete Toda system on the skew square lattice are thus obtained 
from those of SectEl by setting (3 = 0. The discrete Lagrange function obtained by 
summing up the elementary Lagrangians along all edges of a time slice, is equal to 

C(X,X) = 22(A(x k , x k ;e) + A(x k , x k+1 ; -e)) . 

k 

The equations of the Toda system read: 

f(x, x; e) + f{x, z_i; -e) + f{x, x; e) + f(x, x X ] -e) = 0. (49) 

The Hamiltonian form of these equations is 

P=-f(x,x\e) + f(x,x- 1 ;e), p = f(x,x;e) - f(x,Xi,e). (50) 

Proposition 12. The system UHty admits a zero curvature representation with the 
matrices 

L = L(x,p; A) = e p M(x, 0; A) + e~ p M(x, 0; -A), (51) 

V = V(x, e; A) = U(u, u_x; A, A - e). (52) 
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Rewriting (}4*9*j) in the multiplicative form, 

F(x,Xi,e) F(x,x;e) 
F(x,x;e) F(x, x_i;e)' 

and then performing transformations like those mentioned at the end of Sect. |HJ we put 
the discrete Toda system in the rational form: 

U\ — u u — U_i 

In this form it was given in 0]. 

11 Elliptic Toda lattice 

Performing the limit e — > in eqs. (JHUJl (or (3 — > in eqs. (jHEJ), (jHZ| ), we find: 

P = log^^, (53) 
x — 1 

p = -C(s + xi) + C(xi - x) - ((x + x-i) - C(x - + 2((2x). (54) 
From ()5l-}jl . (jHlj) there follow the Newtonian equations of motion: 

x 

— — - = ((xi + x) - C(xi - x) + + x_i) + C(x - - 2£(2x). 
x^ — 1 

This is the elliptic Toda lattice as given in |12j . The rational form of this equation (with 
u = p(x), as usual) was discovered earlier in [221 12E] and reads: 

u-r'{u)/2 _ 1 1 
u 2 — r(u) u — u\ u — U-x 

The matrix L the zero curvature representation of Proposition IT2l does not depend on 
e. As a consequence, this matrix serves also for the continuous time system. 

Proposition 13. The system \5§) . \54\ l admits a zero curvature representation \3ty) with 
the matrices L from / T57|) and A = A(u, U-±; A) = dU(u, A, A — s)/de\ £= o. 

An alternative zero curvature representation was given in J2] (with a less local matrix 
L depending on x,p and Xi). The relation of both representation remains unclear. 

12 Square lattice, explicit scheme 

Last, we consider the square lattice with a slicing as on Fig. |HJ obtained from the trian- 
gular one by erasing the diagonal edges. This can be achieved by setting a 2 = a 3 . For 
our particular setting plj) this means that all parameters of the Toda lattice are equal: 
(3k = £■ The parameters of the corresponding quad-graph (finer square lattice shown by 
dashes on Fig. H3J) are a\ = A, a 2 = A — e. 
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X_i X X\ 

Fig. 14: Fields and wave functions on a square lattice 

Correspondingly, the results on the discrete Toda system on the square lattice are 
obtained from the results of Sect. El by putting (3 = e. The discrete Lagrange function is 
equal to 

£(X,X) = ^2(A(x k , x k ;e) + A(x k , x k+1 ; -e)). 

k 

The equations of the discrete Toda system read: 

f(x, xi, -e) + f(x, x; e) + f(x, -e) + f(x, x; e) = 0, (55) 

and their Hamiltonian form is 

p = -f(x,x;e) + f(x,x 1 ;e)+f(x,x-i;e), p = f(x,x;e). (56) 

The rational form of (|53|) reads: 

p(u\,u, u;e) p(u,u,U-i,£) 
U\ — u u — U-i 

The zero curvature representation is of the form (|29j) . with 

L = L(x,p, A) = e p M(x,e;X) +e- p M(x,-e;-X) 1 
V = V(x, A) = U (u, U-i] A, A — e). 

Notice also that, in virtue of ([S^)l. p = f(x, x; e), hence L = U(u, u; A, A — e). 

It is important to mention that both maps (joTIj) and (JBTIj) approximate in the continuous 
limit e — > the elliptic Toda lattice (|5Hjl . (jHlj) . The first (implicit) map is a Backhand 
transformation of ETL itself (it shares the integrals of motion and commutes with ETL). 
However, the second (explicit) map is a Backlund transformation of ERTL with the 
parameters ft = e. 
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A List Q of integrable equations on quad-graphs 



The first three items of the list IQ1HQ4I are 

a(uo — u 3 )(ui — u 2 ) — b(uo — ui)(u 3 — u 2 ) + S 2 ab(a — b) = 0, 
a(u - u 3 )(ui - u 2 ) - b(u - u x ){u 3 - u 2 ) 

+ ab(a — b)(uo + «i + u 2 + u 3 ) — ab(a — b)(a 2 — ab + b 2 ) = 0, 



(Qi) 
(Q 2 ) 



b(a 2 - 1)(m Mi + u 2 u 3 ) - a(b 2 - l)(u u 3 + Uiu 2 ) 

+ {b 2 - a 2 )(u u 2 + u x u 3 ) - S 2 (a 2 - b 2 ){a 2 - \){b 2 - l)/(4o6) = 0. (Q 3 ) 

The additive three-leg forms are collected in the table below. (In all cases except for 
(Qi !5=o it is actually more convenient to deal with the multiplicative three-leg form (J2J).) 



f(x 0) xi\a) 



U = 0[X) 



a = p(a) 



(EI 



(EI 



5=0 



(5=1 



<5=0 



(Q: 



(H 



.5=1 



a 



2 x 



x - x x 
Xq — Xi + a 



Xi — a 

1 , sinhfxo — xi + a) 

o log ^T7 T 

2 smh(£ — Xi — a) 

1 lo ( x ° + Xl + a )( :ro ~ Xl + a ) 

2 ^ (xo + xi — a)(xo — xi — a) 

1 sinhf^o + Xi + a) sinh(x — x± + a) 

— log r 

2 sinh(xo + x\ — a) sinh(xo — X\ — oc) 

1 j cr(x + Xi + a)cr(x — X\ + a) 

2 cr(x + X\ — a)<j(xo — x x — a) 



x 



X 



exp 2x 



x 



cosh 2x 



a 



a 



exp 2a 



a 



exp 2a 
p{a) 



According to the classification routine developed in the main characteristic of 
an integrable equation (JTJ) is a certain polynomial r{u). For the equations IQ1HQ4I this 
polynomial is equal to 

1, u, u 2 -5 2 , Au 3 -g 2 u-g 3 , 

correspondingly. Using this information, one easily finds the following sequence of degen- 
erations: 

(^31(5=1 — ► 



I 

m 



3]l<5=0 
I 

ra«5=i 



(OH) (5=o 
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where 

¥1 



(33«=i : 92 



(§3*=i 



— > QQ3D<?=o 

(HI])5=i : 

( ^3 ) 5=0 — >■ ( Q7 ) 5=1 



(Qil5=i -> (Qil* 



M - 

#2 
It - 

U - 



u 
a 



4 

3' 



+ 3' 



4a 1 
(a-1) 2 + 3' 



5-3 = 0, u ->!/«, a 



A 



2eV 



a 



1 + 2ea, e -»• 0, 



u/e, 



1 



ea, u — > — h ew, e 
4 



1 + 2ew, a 
ea, e — > 0. 



2ea, 



8a(a + 1) 
" ~ (a - 1)3 • 
-2a' 3 , 



B Some formulas with elliptic functions 

Some of the most useful formulas for the Weierstrass elliptic functions are: 

a(x + a)a(x — a)a(f3 + j)cr(P — 7) = cr(x + [3)a{x — f3)a(a + j)a(a — 7) 

— a(x + 7)cr(x — 7)cr(a + /3)a(a — (3), 



t(x) + C{y) + ({*) - <(x + y + z) 



1 p(x) p'{x) 

1 p{y) p(y) 

1 p(z) p'{z) 



p{x) - p(y) 



a(x + y)cr(y + z)a{z + x) 
a(x)a(y)a(z)a(x + y + z) 

a(x + y + z)a(x — y)cr{y — z)a(z — x) 
a 3 (a;)cr 3 (?/)o" 3 (2;) 

a(x + y)a(x - y) 



a 2 (x)a 2 (y) 

A fundamental object related to the elliptic curve £ is the polynomial 

H(u, v, w) = (uv + vw + wu + (?2/4) 2 — (u + v + w)(Auvw — #3). 



(B.l) 
(B.2) 

(B.3) 
(B.4) 

(B.5) 



Indeed, the equation H(p(x), p(y), p(x ± y)) — is the addition theorem for the go- 
function. This follows from the formula 



H(p(x),p{y),p(z)) 



<y(x + y + z)a(—x + y + z)a(x — y + z)cr(x + y — z) 
a\x)a\y)a\z) 



(B.6) 



There holds the identity H 2 — 2HH VV = r(u)r(w). This gives rise to a one-parameter 
family of symmetric biquadratics h(u,v; (3) with the property h 2 — 2hh vv = r(u): 



h(u, v- (3) = H(u, v, b)/ v^fej = H(u, v, p((3))/p'(P). 



(B.7) 



Recall H that the polynomial h is related to the basic polynomial Q from eq. (0) by the 
formula Q U2 Qu 3 - QQu 2 u 3 = k(a, (3)h(u , u\\ a). 
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To an arbitrary quadratic polynomial h(u) = c 2 w 2 + 2c\U + c there corresponds the 
symmetric affine- linear polynomial p(u,w) with the property h(u) = p(u,u), namely 
p = C2UW + C\(u + w) + Co- In a more invariant fashion, 2p = 2h — (u — w)h u . For the 
polynomial h = h(u,v; (3) from ()B.7jl this polynomial p is given by: 

p(u,v,w;[3) = {M ~^ ))2 (uw-(p(y+P)+p(y-P))^^+p(y+P)p(y-P)), (B.8) 

where, as usual, u = p(x), v = p(y) and w = p{z). This formula follows easily from 

h(u, v- (3) = tefo)-P(fl) 2 ( M _ p ( y + /3))(„ - p(y - fl), 



which, in turn, is a direct consequence of (jBJy), f)B.7j) . 
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